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We explore the crucial role of relative space-time positioning between the two detectors in an 
operational two-party entanglement-harvesting protocol. Specifically we show that the protocol 
is robust if imprecision in spatial positioning and clock synchronization are much smaller than 
the spatial separation between the detectors and its light-crossing time thereof. This in principle 
guarantees robustness if the imprecision is comparable to a few times the size of the detectors, which 
suggests entanglement harvesting could be explored for tabletop experiments. On the other hand, 
keeping the effects of this imprecision under control would be demanding on astronomical scales. 


I. INTRODUCTION 

The vacuum state of a long-range field is unentangled 
with respect to its extended nonlocal modes but can be 
entangled with respect to local modes [I] corresponding 
to localized detectors. Extracting this naturally endowed 
resource from a quantum field vacuum via localized de¬ 
tectors, such as atomic probes, is called vacuum entan¬ 
glement harvesting EHU- The process can be understood 
in the following terms: An actuator turns the interaction 
between these two detectors and the field on and off, 
resulting in two initially unentangled detectors evolving 
into an entangled pair, after ignoring, or tracing out, the 
field itself. This process is valid even for two spacelike- 
separated detectors. 

The resultant pair of entangled detectors can in prin¬ 
ciple serve as a resource for performing quantum infor¬ 
mation tasks such as teleportation 00, superdense cod¬ 
ing muni, or fingerprinting [uina. This resource is es¬ 
pecially useful when enhanced by entanglement-farming 
protocols using successive pairs of detectors |T3]. We 
are particularly interested in the role and limitations to 
this protocol that arise form the requirement to control 
the spatial separation between the detectors and their 
switching synchronization. 

Previous studies do not consider the spacetime posi¬ 
tion problem but rather assume that separation is known 
and a precise clock is shared between the two parties 
(henceforth Alice or A and Bob or B) undertaking the 
entanglement-harvesting protocol. Here we establish the 
sensitivity of entanglement harvesting under imprecision 
in spacetime positioning and determine when the proto¬ 
col is robust under such imprecision. 

The usual model for entanglement harvesting treats 
two localized Unruh-DeWitt detectors [U [14] linearly 
coupled to a scalar field in a vacuum state according to 
well behaved switching functions nans]. Despite its sim¬ 


plicity, the Unruh-DeWitt model successfully captures 
essential features of the light-matter interaction mm- 
We treat the field as occupying a (3 + 1 (-dimensional flat 
geometry with spacetime coordinate x := ( x,t ). Bipar¬ 
tite entanglement can be assessed by calculating negativ¬ 
ity Em. a positive value of negativity is a necessary and 
sufficient condition for entanglement in the case of two 
qubits [ZCTt I2T1 . 

II. ENTANGLEMENT HARVESTING 
PROTOCOL 

An Unruh-DeWitt detector is a localized particle cen¬ 
tered at x with spatial smearing function F{x) (e.g., 
spatial wavefunction of electron for single-atom detector) 
and has two distinct internal energy levels | g) and |e) for 
ground and excited state, respectively, separated by an 
energy gap D. For |e) = <j + | g), working in the interac¬ 
tion picture, the detector’s monopole moment is [17j [IS] 

fi{t) =a + e int + a-e~ int . (1) 

The detector monopole’s moment is linearly coupled to 
a massless scalar field, whose operator in the interaction 
picture is given by 

(2) 

for k := (fc,w = |fc|) with H.c. designating Hermitian 
conjugate. We denote the fc th -mode annihilation opera¬ 
tor by a(k ), and its adjoint creates the plane-wave single¬ 
photon Fock state 

ot(fc)|0> = |l(fc)> (3) 

from the vacuum |0). 
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Given field-detector coupling strength A and switching 
function x(t), the interaction Hamiltonian is 


Hi = Xx{t)p,(t) / d 3 xF(x)<j>(x, t ), 


( 4 ) 


where F(x) is the spatial profile of the detector. For 
a non-pointlike profile, the behaviour of the detector 
is regular even for piecewise continuous switching func¬ 
tions [HI. We introduce the detector’s form factor 


F v (k)= <Tx(2n )-iF„{x)e lkx , 


( 5 ) 


where v £ {A,B} labels Alice’s and Bob’s detectors. 
From this, we write 

Hi= ^2 Kxv{t)iiv{t) 

v=A,B 

x / (k> h{k> +H,c -) ■ (6) 
and the perturbative evolution operator is 

/ OO /*00 nt 

dtHi(t)- d t dt'Hi{t)Hi(t')+0(X 3 ) 

-oo J —OO J —OO 

s 

( 7 ) 


um 


( 7 ( 2 ) 


for 


A := max Ay <!!:= min fly. 


( 8 ) 


Inter-detector separation and detector switching can 
be partitioned into two regimes: lightlike contact, 
whereby real quanta can be exchanged, and spacelike sep¬ 
aration, whereby they cannot be exchanged. Notice that 
it is possible that the two detectors communicate when 
they are in timelike contact even though energy does not 
necessarily flow from Alice to Bob [22H2o] . In the space- 
time diagram in Fig. [I] we present the possible cases of 
relative positioning of Alice’s and Bob’s detectors, fol¬ 
lowing the same conventions as in [231 EH] ■ 

The final state for the detectors and field is 

p = Po + p (1) + p (2) + 0(A 3 ) (9) 

for initial ground state 

Po ■= |0 gg) {Ogg\ (10) 

and 

P {1) =U {1) p 0 + p 0 U (1) \ p (2) = p[ 2) + p2\ 
p {2) :=U^p 0 U^\ p< 2) := U^po + p 0 U^ . (11) 

The final reduced two-detector state, 

Pab := trf p (12) 
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Figure 1. (Color online) Spacetime diagram showing the pos¬ 
sible relative placement of the switching periods of Alice’s and 
Bob’s detectors. We follow the same conventions as in 1231 1261. 
showing the following possible relative position of the detec¬ 
tors Alice and Bob: 1-purely spacelike, 2,3,4-(partially) light 
connected, 5-purely timelike, 6-lightlike only for part of Bob’s 
interaction. 


for trf denoting partial trace over the field, satisfies 

ti'f (t/ (1) po) = 0 (13) 

so the leading-order change to pab is 0(A 2 ). 

Let Alice and Bob employ identical detectors: F{x) 
and A are the same. We assume a real even spatial profile 
F{x) = F(—x), so 


F(k) = F(—k). 

We introduce shifted frequencies 

Wy (ft) := w(fc) ± f \ 

and integrals 


(14) 

(15) 


/y(fc) := Ay 


i k-x u 


\phjj 


F(k) / d tXuMertW (16) 


and, for r = x a — (and r := |r|), 

/ OO nt 

d t / d t' Xa^Oxb (t) 

-oo J —oo L 


x e ' 




m 2 > k .r 

2co 

(17) 


Then we define 


Xyy' := j d 3 ki:(k)I„,(k), J ■= Jd 3 kJ(k), (18) 


and 


2^ := Taa ±Ibb, 


(19) 


and these integrals (18) suffice to characterize pab fully. 
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Straightforward but lengthy algebra leads to 

p W= J d 3 k J d 3 k'(l A (k)I A (k')*\l(k)eg)(l(k')eg\ 

+ I B (k)I B (k')* \l(k)ge)(l(k')ge\ 

+ I A (k)I B (k')* \l(k)eg)(l(k')ge\ 

+ I B (k)I A (k')* \l(k)ge)(l(k')eg\ ) (20) 

with double excitations achieved via counter-rotating 
term a + aJ(k ), and higher states are not excited to this 

/ o\ 

perturbation order. A similar calculation yields p\ . To 
compute pab, we trace out the field to obtain 


tr fPi 2) = / d 3 fc | / A(fc)| / \eg)(eg\ + |/ B (fc)| 2 \ge){ge\ 


+ (l A (k)IUk / )\eg)(ge\+K.c.) 


trfp^ = — / d d fe 


2ReC(k) \gg)(gg\ 


+ J( k ) \ ee )(gg\ + J*i k ) \gg)( ee \ 


( 21 ) 


and / d 3 fc ReC(fe) conveniently does not need to be cal¬ 
culated as it can be obtained from the other terms as each 
perturbative correction is traceless order-by-order m- 


Pab = 


{| gg) ■ 

. I ge) , 

\ e g) 1 

|ee)}, 

-i+ 

0 

0 

-J*’ 

0 

Z BB 

Iab 

0 

0 

7-* 

X AB 

Iaa 

0 

-J 

0 

0 

0 


+ 0(A 4 ). (22) 


In contrast to one of the original entanglement-harvesting 
proposals by Reznik et al. 2], one diagonal term van¬ 
ishes here because exciting both detectors requires 0(A 4 ) 
terms. If we were to include this term, we would also have 
to include all the additional 0(A 4 ) terms and not only the 
double-excitation one. (Notice, however, that as Reznik 
et al. calculate entanglement to second order, all their 
entanglement results remain correct in spite of this small 
inconsistency in the Dyson expansion.) 


III. BIPARTITE ENTANGLEMENT 


For a two-qubit system, negativity N conveniently es¬ 
timates bipartite entanglement [T9]. For transpose oper¬ 
ator T, application of the partial transpose operation to 
a general density matrix yields 


1®T : Yl p ki I*) 0'M fc ) (1 ^ I*) 01® 10 ( fc i ( 2s ) 

ijkl ijkl 


whose negative eigenvalues sum to A f. 


From (22 I, 


(1 <8> T) pas 


/1-J+ 0 

[ 0 I BB 

0 -J 

V x* AB 0 


0 Zab\ 
-J* 0 

Taa 0 I 

0 0 / 


+ 0(A 4 ) 
(24) 


with only one negative eigenvalue at second order: 


M = ~, 


X+ -^(T-)2 + 4 |Tr 


+ 0(A 4 ). (25) 


For identical detectors switch ed on for the same amount 
of time ( X~ = 0), negativity (25) simplifies to 

N=\J\-1. (26) 

As pointed out by Reznik et al. [3j, intuitively the non¬ 
local term J oc X A X B must exceed local noise I + oc A 2 
for entanglement to arise. 

We consider the case of sudden switching; concomitant 
ultraviolet divergence is offset by spatial smearing [2T 
( 281 . Switching is described by the rectangular function 


1, if |t| < 1/2, 
n (t) = { 1 / 2 , if \t\ = 1 / 2 , 

0, if \t\ > 1/2 


with Fourier transform 
/>00 

/ dt n (t)e~ iut 


= si “ (I) 


The detector is switched on at T and off at T' so 

't-T 1' 


x(t) = n 


for 


T± := 


2 T- 


V ±T 


(27) 


(28) 


(29) 


(30) 


As 


piojT' _ p i ujT i 

____ =- e k 

ui 


I T‘ 


LO 


=2i T e luT sincwT , 


(31) 


we obtain 


Q i k-Xjy 


I v (k) =2X V T~ e iu ”W T " —j= 
V2 co 

x sinc(oj j (k)T~)F(k) 


(32) 


so 


T 

J -W 


d 3 fc sin 2 
2 ui 


^+(ft)r; 


( 2F(fc) 

yw1"(fc) 


2 


(33) 


Making use of the compact support of y„(t) over [T v , T'], 
we can write 

J(k) = X A X B e ikr (j AB (fc) + J BA (fc)) (34) 

where 

r T l r l 

J vv ’(k) := / dt dt' Xu ,(t')e^' {k)t '-^ (k)t) (35) 

J -oo 

with v, v' € {A, B}. 

The three distinct timing regimes are 
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(i) T' < T v > with (k) = 0; 

(ii) T', < T v for which 

3 vv >(k) = - 4T-T ! ^e i K' (fc)T ^- w - T - +(fc) ) 

x sine co~(k)T~ sine wJ(fe)rJj; (36) 

and 

(iii) the detector is on for overlapping times for which 
we assume without loss of generality that T„/ < T v 
and T' v , < T' v . 



r (x lOOcr) 


As ,J v ,y ( k ) ^ 0 only if a detector is on, we consider only 
the no-overlap time domain 


T 0 := [T„i, T v ) U {T ' v ,, T' v \ 
or the overlapping interval 

Ti := [T V ,T' V ,]. 

Thus, 

JVII 1 Ayy' -A ui 

with 

t ' i T 
" ( k)t v 


(37) 

(38) 

(39) 


it o u ( k)t 


Kuv{k) =- 


u>„{k) w+,(fc) 

T' T > 

iui+ ( k)t 
T', 


a —i uj u ( k)t 


over 7o and 


x*Ak) =- 


u v {k) w+,(fc) 

i(u+(fc)T„-wy(fe)T',) 


(40) 


(fe)wj,(fc) 

u;+(fc)e ir - n ^-u,j(fe)e iT ^ n ^ 

w^(fe)o;+,(fe)n+ I/ , 


over 71 for 


:— fly + fly'. 


(41) 


(42) 


We specialize to Gaussian spatial smearing of the de¬ 
tectors, 


F(x) = 


exp (~\x\ 2 /a 2 ) 

(^) 3 


(43) 


which implies negligible overlap between smeared de¬ 
tectors for reasonable separations r a. For non¬ 
overlapping switching functions, Eqs. (33)-(36) and (18) 
yield 


X = £ 

- L 'VV — 9 

7 T z 


dw 


1 , 2 2 
uie * 


Wi 


+2 


sin^ I '-fT- 


\J\ = 


AaAb 

47r 2 r 


^ —t 


_ 1 2 2 
dcusincere !UIT ■ 


T' 

R 


B . + , 

e lw A 4 


(44) 

Ti 

T a 


Wf 




(a) 



(b) 

Figure 2. (Color online) (a) Correlations term \J\ and 
(b) projection of pab onto |$ + ) vs distance r with fl = 1 
for detectors with 2T^" B = lOOcr, Tb — T' A = 50<r, and 
<7 = 0.001, Vertical (red) dashed lines indicate the start and 
end of lightlike contact for duration 10er. Both peak on light- 
contact showing that the maximum of correlation between the 
atoms occurs at light-contact. Notice the leakage outside of 
the light-cone. 


with important quantities 

u;a = Ifcler 


(45) 


being the wavenumber-spread product and uir being the 
normalized separation in terms of wavenumber. As J is 
asymptotically proportional to r^ 1 , entanglement, which 
involves competition between r-independent T vl , vs J7, 
decreases as r increases. 

Although for simplicity and clarity we have chosen for 
our study to employ a sudden switching function, no¬ 
tice that the harvested entanglement in the case of a 
sudden switching is noticeably smaller than in the case 
of Gaussian switching |6]. Note, however, that for any 
kind of switching, the fundamental feature which leads 
to entanglement harvesting (at leading order) is that the 
non-local term \J\ dominates over the local noise X in 

Eq. (1261). 


IV. SENSITIVITY TO SPACE-TIME 
POSITIONING 

As is well understood (See, e.g. [6]), harvested entan¬ 
glement is evidently maximized for the two detectors be- 
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ing lightlike-separated but can be nonzero even for detec¬ 
tors which remain spacelike separated. One can clearly 
see this by looking at the non-local term \J |, which peaks 
if the two detectors are on for a long time while light- 
connected but falls with increased r. This decline with 
respect to r applies to spacelike entanglement harvesting 
for which detectors do not exchange real photons and 
the correlation term \ J\ (and, therefore, the harvested 
entanglement) deteriorates due to decay of vacuum cor¬ 
relations with r. 

Additionally, for any entanglement in pab to be useful, 
Alice and Bob need an accurate description of the mu¬ 
tually shared (mixed) state pab, which can be expressed 
in the Bell basis 


|$±) =(l,0,0,±l) T / v / 2, 

|d/±) =(0,l,±l,0) T /v / 2. (46) 

The state’s Bell-state fraction for the final state of the 
detectors is 

(4> ± |pAB|$ ± )=i-i[l-T2Re( l 7)], 

(4^1 Pab |^ ± > =1 [ 1 + ± 2Re(X AB )] ■ (47) 

In Fig. [2] we show both the behaviour of the correlation 
term \J\ and the projection of the shared correlated state 
onto | c f> + ). Evidently viable entanglement harvesting is 
maximized by the inter-detector light-crossing time being 
order unity expressed in terms of r in natural units. 

To exploit harvested entanglement, Alice and Bob re¬ 
quire an accurate description of the r-dependent shared 
density matrix pab, but r necessarily has nonzero un¬ 
certainty 8. We treat the joint state, subject to a mean 
separation r 0 with uncertainty of 8, as a smeared state 

/' e -(r-r 0 ) 2 /S 2 

Pab{$) = J drPr(r)p AB (r), Pr(r) = --, (48) 

with the Gaussian distribution Pr(r) formally permitting 
an unplrysical negative distance, which is not problematic 
for reasonable separations. 

Importantly, in the matrix elements of pab (8 ), the only 
deleterious effect on A f arises in the term J as the terms 
are r-independent as expected from their local na¬ 
ture. This difference allows us to focus the study of how 
synchronization affects the ability of Alice and Bob to 
harvest entanglement exclusively on the behaviour of the 
correlation term \ J\. 

For the simpler non-overlapping case of the switching 
functions, we can calculate a closed expression for \J\ 
modulo the integral over field frequencies using the fact 
that 


dr sin(|fc|r) 


(r-r n ) 2 


~(ro\ 2 r 
= 7re t 6 > 1m 


Erfi 


ro 

8 


i(5|fc| 

2 


( 49 ) 


where Erfi(z) = —iErf(iz) is the imaginary error func¬ 
tion. Taking expression ([49]) into account, the effect of a 
relative-positioning error 8 on the correlation term \J\ is 


\J{5)\ =^^ e -Po/5) s 


4(5 7T 3 



„ „ (r n i(5aA 

, X 

/ dcj Im 

Erf i -r + ir) 

J 0 

L W 2 )\ 


t 


1 / , 2 r- r 2 

x e“ 5 - 


T' 


B . + , 


t b 


Tf 


Ta 


Wf 


W 


(50) 


We see from Eq. (501, that, if the uncertainty of the 


inter-detector separation is fixed, | J\ decreases with de¬ 
creasing spatial separation ro between the atoms. Specif¬ 
ically, the Gaussian term in Eq. ( |50| ) can dramatically 
reduce entanglement-harvesting capability for long dis¬ 
tances. On the other hand, if 8 increases, we see that in 
the limit that 8 is larger than all the other length scales in 
the problem (5 —► 00 ), both the Gaussian factor and the 
imaginary part of the Erfi function in Eq. (50) asymptot¬ 
ically approach unity, and then \J\ ~ (5 -T_ ^ This shows 
that entanglement harvesting is impacted by spatial un¬ 
certainty in the positioning of the detector, and the decay 
is just of a polynomial kind. 

If the two detectors held by Alice and Bob are fixed 
precisely, within precision cr for each detector, prior 
to commencing the entanglement-harvesting protocol, 
and Alice and Bob are informed of the value ro a to a 
high degree of precision, then the effect of nonzero 8 can 
be ignored. If, on the other hand, Alice and Bob need to 
determine their separation as the first stage of the pro¬ 
tocol, which is a reasonable requirement for harvesting 
entanglement from the cosmological vacuum, then estab¬ 
lishing separation to precision within 0(a) independent 
of r is daunting. 

We now analyze the resultant uncertainty due to fun¬ 
damental imprecision in synchronizing their respective 
reference frames. Due to the covariant nature of the 
light-matter interaction, the impacts on the atom-field 
dynamics of a distance uncertainty and of a relative time 
uncertainty are the same (as one of them becomes the 
other under an appropriate reference frame change). This 
rule is true if the time uncertainty is insufficiently small 
to resolve the additional time scale in the problem that 
is not related, in a Lorentz-covariant way, to a physical 
length scale (namely, S2 _1 ). The reason the scale fl -1 is 
different from the others is that the internal dynamics of 
the atoms is considered (to a very good approximation) 
to be non-relativistic in the Unruh-Dewitt model (and 
any other usual models of light-matter interaction). 

Although a lengthy task, numerically checking verifies 
that, for time uncertainties exceeding f2 _1 , the uncertain¬ 
ties in space and time impact on entanglement harvesting 
in the same way as intimated by the above reasoning. 
Thus, we focus our detailed analysis on the impact of 
a spatial uncertainty 8 on harvesting, which we also call 
‘synchronization error’ for uncertainties much larger than 
the timescale fl -1 . As a last comment, the scale of fl -1 
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Figure 3. (Color online) Correlation term \J\ (responsible 
for entanglement harvesting) vs separation imprecision 8 for 
ro = 150 ct, 2T^ B = 100 ( 7 , Tb — T' A = 50a and a = 0.001 
evaluated at optimal separation and switching time accord¬ 
ing to Fig. [2ja) . The vertical dashed (red) line represents 
the distance between the two detectors. Although the phys¬ 
ically relevant region covers only the left of the red dashed 
line (uncertainties in positioning below the separation of the 
atoms, 5 < r) we plot the large uncertainty regime to show 
the polynomial decay of \J\ with the imprecision 8. 


for the Hydrogen atom’s first transition (of wavelength 
A ~ 100 nm) is fl _1 = c _1 A ~ 10 -16 s. The assump¬ 
tion that the main sources of uncertainty are going to 
be larger than this timescale is, therefore, a good ap¬ 
proximation, and henceforth we consider that imprecise 
determination of r and synchronization limitations are 
alike. 

Alice’s and Bob’s harvestable entanglement from the 
field vacuum is robust under imprecision as exemplified 
in Fig. [3] where we plot the relative decrease in the cor¬ 
relation term as the imprecision in the determination of 
r is increased. In particular we plot the ratio 


m 


\m\ 


(51) 


evaluated in the vicinity of the peak in Fig. [2j We see 
that the correlation term |J7’(^)| decays 60% when the 
uncertainty in the determination of relative distance is 
the separation of the detectors. Notice that in most ex¬ 
periments the uncertainty in the separation between the 
two atoms would be smaller than the distance between 
the atoms r, and typically larger than the atomic size er, 
that is r > 8 a. This means that the most physically 
relevant region of the plot is to the left of the red dashed 
line. We plot the behaviour of the correlation term |d%(<5) | 
for uncertainties larger than the separation between the 
detectors to illustrate the slow (polynomial) decay of the 
ability of the setup to harvest entanglement. 

When the uncertainty 8 is larger than the distance be¬ 
tween the detectors, one may think that Fig. [2] may con¬ 
tain non-physical situations where the two detectors are 


overlapping. However, as r <C <r, if we think of the uncer¬ 
tainty emerging from an ensemble of experiments where 
the positions of the centre of mass of the two detectors 
in every experiment of the ensemble are randomly set 
along the Gaussian-distributed length 8, the probability 
of the two detectors landing in a region where they would 
have a non-negligible overlap (of the order of a few a) is 
exponentially suppressed. Therefore, the effect of these 
possible superpositions of the two detectors on features 
of the plot in Fig. [2] is negligible. 

It is known that, even in the cases where entangle¬ 
ment harvesting is possible for a broad set of physical 
parameters (See, e.g., [S]), the correlation term \J\ never 
becomes much larger than the noise term X. Instead, in 
all scenarios for which entanglement harvesting is possi¬ 
ble, the correlation term becomes larger but of the same 
order of magnitude than the noise term. 

For entanglement harvesting, \ J\ > X is crucial, and 
damping by 50% of the correlation term already suf¬ 
fices to prevent entanglement harvesting from happening. 
Thus, according to our results, if the uncertainty in the 
determination of the spacetime synchronization is of the 
order of the separation between the detectors, entangle¬ 
ment harvesting is impeded. 

As we can see in (50), the decay of the term \J\ is 
dictated by the ratio of the spatial separation to the 
distance uncertainty (ro /8) and not the smearing length 
scale of the detectors a. This means that, even though 
in Fig. [3] we show regimes where 8 cr, the decay in 
the non-local term \J\ is entirely due to the increase of 
uncertainty 8 and not to any kind of overlap between the 
detector smearings or any other effect related to the size 
of the detectors relative to the magnitude of their spatial 
separation. 


V. CONCLUSIONS 

We have reviewed that entanglement harvesting is the 
result of the the non-local terms \J\ dominating over the 
local noise terms X in the interaction of a system of two 
detectors with the vacuum state of a quantum field. We 
have analyzed how imprecisions on the space-time syn¬ 
chronization of the two detectors (their relative distance 
and their ability to synchronize their clocks) affects both 
the local noise and the non-local terms. 

We have shown that while, of course, the local noise 
terms are unaffected by imprecision in the synchroniza¬ 
tion of the detectors, the non-local term gets exponen¬ 
tially damped in the square of the relative uncertainty in 
the space-time positioning. This means that in scenarios 
where entanglement harvesting is possible, the harvesting 
of entanglement between two particle detectors is robust 
under relative imprecisions in the spatial and temporal 
positioning if they are below the scale of their spatial 
separation and the light crossing time between them. In 
contrast, the ability to harvest entanglement deteriorates 
rapidly when the precision in the space-time positioning 
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becomes non-negligible with respect to the separation be¬ 
tween detectors. 

If the precision in spacetime positioning is comparable 
to the size of the detectors a <C r 0 , as could be the case 
of atomic-scale detectors in a well controlled table-top 
experiment, our results hint that vacuum entanglement 
harvesting would be robust under distance-time measure¬ 
ment imprecisions. On astronomical scales, however, in 
situ measurements of distance are generally not possible 
so distances are instead inferred from various measured 
quantities such as the product of redshift and Hubble’s 
constant [2S] where even establishing the value of the con¬ 
stant is challenging m- This limitation indicates that 
our results, when extrapolated to a cosmological scenario 


where the field is conformally coupled to curvature, could 
be particularly relevant when considering entanglement 
harvesting in cosmological contexts |23| EB. 
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